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We consider an admissible Riemannian polyhedron with piece-wise smooth 
boundary. The associated Laplace defines the boundary spectral data as the set of 
eigenvalues and restrictions to the boundary of the corresponding eigenfunctions. 
In this paper we prove that the boundary spectral data prescribed on an open sub- 
set of the polyhedron boundary determine the admissible Riemannian polyhedron 
uniquely. 

1 Introduction 

Recent years have seen some very significant achievements in the study of 
inverse boundary-value problems in a single component body. Mathemati- 
cally such body is described by a PDE or a system of PDE's with relatively 
smooth coefficients. Starting from the pioneering works [7J and [42j, inverse 
boundary-value problems were solved, at least on the level of uniqueness and 
sometimes conditional stability, for a wide range of scalar inverse problems, 
both isotropic and anisotropic, see e.g. @], [8], [S], [20], [23], |2§], [30], [21], 
[32] , [35] . [38], [JT] for a far from complete list of references, with further ref- 
erences in monographs [17J or [2T]. Moreover, for such media there appeared 
a number of important results in the study of the inverse boundary-value 
problems for systems of PDE's corresponding to physically important mod- 
els of electromagnetism, elasticity and Dirac equations, see e.g. [18], [27] . 
[28], [33], [M], [36], [37J. 

Much less is known, however, about the inverse boundary-value problems 
for a multicomponent medium. Mathematically, such medium is described 
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by PDE or system of PDE's with piece-wise smooth coefficients with dif- 
ferent subdomains of the regularity of coefficients corresponding to different 
components of the medium. Clearly, the study of inverse problems for the 
multicomponent media is of substantial importance for practical applications. 
Imagine, for example, a human body with bones, muscle tissue, lungs, etc. 
each of those having distinctive values of material parameters, or an upper 
crust of the Earth which is a composition of clay, sand, rock, oil, water, etc. 
A complete answer to the inverse boundary problems in a multicomponent 
medium, at least when the data are measured on the whole boundary, is ob- 
tained only for the two-dimensional case. Namely, it was shown in [2], [3] that 
the Calderon inverse boundary problem in the 2D case has a unique solution 
in the class of L°°— coefficients. Clearly, these results cover also the case of 
a multicomponent medium. In higher dimensions, the results are restricted 
mainly to the inverse obstacle problem. In these problems the goal is to find 
a shape of an inclusion inside a given medium which parameters are known 
a priori. In the case when parameters of a medium and/or inclusion are 
unknown they are assumed to be homogeneous throughout each component, 
see e.g. [TJ, [H], [T5], [T5], [20] ■ Having said so, we should note that there 
exist powerful methods to find singularities for coefficients of lower order, see 
e.g. [I!]. 

This paper is devoted to the study of the inverse boundary spectral prob- 
lem for the Laplace operator in a multicomponent medium. To be more 
precise, we assume that the domain occupied by the medium consists of 
a finite number of subdomains with piece-wise smooth boundaries between 
them. The metric tensor in each subdomain is smooth but does have jump 
singularity across the interfaces, i.e. the boundaries between adjacent sub- 
domains. Adding proper transmission conditions across the interfaces and 
boundary conditions on the domain's boundary, defines a Laplace operator 
which, from the spectral point of view, has effectively the same properties 
as the Laplace operator in a single component medium. Mathematically, the 
considered medium may be described as a Riemannian polyhedron. Leav- 
ing exact definitions of an appropriate Riemannian polyhedron to the next 
section, imagine an n— dimensional simplicial complex Ai where simpleces 
can be glued together, pairwise, along their (n — 1)— dimensional faces which 
we continue to call interfaces (sometimes (n — 1) — interfaces). Imagine now 
that each simplex has its own smooth metric g which, in principle, may have 
jumps across interfaces between adjacent simpleces. This, together with some 
additional geometric/comb inatoric conditions described in section [21 defines 
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a Riemannian polyhedron (A4, g). Starting from the corresponding Dirichlet 
form on H 1 (Ai, g)— functions and using standard methods of spectral the- 
ory, the Laplace operator with Neumann boundary conditions, A, is then 
well-defined in L 2 (Ai,g). Denote by {Xk, <fk}k%i ^ ne se ^ °f a ^ eigenvalues, 
counting multiplicity, and corresponding orthonormal eigenfunctions of A. 
Let T C dM be open. 

Definition 1.1 The collection (T, {Xk, c^fc|r}fcLi) is called the (local) bound- 
ary spectral data (LBSD) of the Riemannian polyhedron (A4,g). 

Let now (Ai , g) and (M. , g) be two Riemannian polyhedra with LBSD 
(r, {A fc , ^fc|r}fcLi) and (f, {X k , ^ fc |f}£Li), correspondingly. 

Definition 1.2 LBSD for (M,g) and (M,g) are equivalent if 

1. T and T are homeomorphic, x : Y — > T; 

Afc = Afc, A; = 1, 2, ... ; 

3. // Xk has multiplicity m + 1, m — 0, 1, . . . ; i.e. Xk = Xk+i = Xk+ m , then 
there is an (m + 1) x (m + 1) unitary matrix Uk such that 

(ip k \r, V^fc+mlr) = U fc (x*if> k \f, ■ ■ ■ , x*<Pk+ m \f)- 

We can now formulate the main result of the paper: 

Theorem 1.3 Let (A4,g) and (Ai,g) be two admissible Riemannian poly- 
hedra. Let, in addition, the metric tensors g and g do have jumps across 
all (n — 1) — interfaces in M. and Ai, correspondingly. Assume that LBSD 

(r, {Afc, <£fc|r}fcli) and (t, {Xk, (fkl^k^ij are equivalent. Then (Ai,g) and 
(Ai,g) are isometric. 

Let us make some comments on this theorem: 

1. If Qi is an n— dimensional simplex of Ai with a smooth metric gi, then 
gi determines an inner metric on any / dimensional, / < n, simplex 
of Ai which lies in Qi (here and later we assume each simplex to be 
close). In particular, any (n — 1)— interface 7 of Ai belongs to two 
adjacent n— simpleces, which we often denote in such case f2_ and Q + , 
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and therefore, has two different metric tensors g_| 7 and g+\ 7 . By a 
metric tensor g having a jump singularity across 7 we mean that, for 
any p <G 7, 

9-(p) ^ 9+(p)- 

This assumption is of a technical nature and, in section [6] we will sig- 
nificantly weaken it. 

2. As shown in section [2j any Riemannian polyhedron has a natural struc- 
ture of a metric space. The isometry of (Ai,g) and (Ai,g) is under- 
stood with respect to these metric structures. 

3. The boundary dA4 of a Riemannian polyhedron (Ai, g) is itself a (n — 
1)— dimensional Riemannian polyhedron, probably disconnected. As 
T, T are open subsets of dM., dA4, by reducing them if necessary we 
assume that T and T are open subsets of some {n — 1)— dimensional 
simplex of dAi, dAi, correspondingly. In the future, we will always 
assume this condition to be true. 

The plan of the paper is as follows: In section [2] we provide some prelim- 
inary material on geometry of Riemannian polyhedra and properties of the 
Laplace operator on them. Section [3] is devoted to the description and some 
properties of the non-stationary Gaussian beams on a Riemannian polyhe- 
dron. We prove Theorem II .31 in sections H] and [51 The last section [6] is devoted 
to some generalizations and open questions. 

2 Preliminary constructions 

2.1 Admissible Riemannian polyhedron 

In this section we will introduce, following mainly [10] and [6], an admissible 
Riemannian polyhedron which is the main object of the paper. We start with 
a closed n— dimensional finite simplicial complex 



2 PRELIMINARY CONSTRUCTIONS 



5 






(A) 



Figure 1: Case (A) is prohibited because its structure is not (n— 1)- chainable; 
Case (B) is prohibited as it is not dimensionally homogeneous; Case (C) is 
appropriate 

where fli are closed n— dimensional simpleces of Ai, with f2* nt standing for 
the interior of Qi which is an open subset of Ai. We assume that Ai is di- 
mensionally homogeneous, i.e. any A;— simplex, < k < n, of Ai is contained 
in at least one We assume also that any (n — 1)— dimensional simplex 
7 belongs either to two different n simpleces, Qi and Qj, which in this case 
we will often denote by f2_ and Q + , or to only one n simplex In the 
former case we call 7 an interface (sometimes (n — 1)— dimensional interface) 
between f2_ and Q + , in the latter case we call 7 a boundary (n — 1)— simplex 
with (n — 1)— simpleces having this property forming the boundary dAi. 
We denote by Ai , < k < n the k— skeleton of Ai which consists of all 
fc— simpleces of Ai. Clearly, Ai = Ai n . We use notations 



Following [10], we assume that Ai is (n — 1)— chainable, i.e. Ai reg is path 
connected, see Fig. [TJ 

Assume now that each n— simplex f2j is equipped with a smooth (up to 
dQi) Riemannian metric i.e. (Qi, gi) is a smooth Riemannian manifold 
with a piecewise smooth boundary. This makes it possible to introduce the 
arclength for admissible paths n : [0, a] — > Ai. We call a path 77 admissible 
if ?]~ 1 (A^ mt ) C [0,a] is a (relatively) open subset of [0, a] of full measure 
and, if rj(a, (3) is in some n— simplex Q mt , then i] : (a, (3) —>■ Q mt is piecewise 
smooth. Naturally, the arclength \r)(a, of the path r\ between 77(a) and 
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T](/3) is taken as 

\ v (a,(3)\= [g m Mt))V m {t)nM'dt, (1) 

J a 

where rjj(t), a < t < (5 are, for example, baricentric coordinates in Q. As 
r]~ 1 (Ai mt ) fl (0, a) consists of at most a countable number of open intervals 

(aci, fli) we define 

|r,[0,a]h^hKA)|. (2) 

i 

Next we introduce, for any p, q G Ai, the distance, <i(p, q), 

d(p, q) = inf \rj\, 
v 

where infenum is taken over all admissible paths connecting p and q. This 
makes (A4, g) into a metric space with its metric topology being the same 
as the topology of a simplicial complex, see [10] . 

Definition 2.1 (A4, g) is an admissible Riemannian polyhedron if, for any 
p,qEM, 

d(p, q) = inf | |, 
v 

where rj run over the subset of admissible paths between p and q such that 
rj' 1 [\jM k ) \({O}U{1})=0. 

\fc=0 / 

As M. is finite, the above condition is independent of a particular choice of 
metric g. 

2.2 Boundary normal and interface coordinates 

In addition to the baricentric coordinates in any Q mt , we will often use bound- 
ary normal or interface coordinates associated with (n — 1)— subsimpleces of 



2 PRELIMINARY CONSTRUCTIONS 



7 




(A) 



(B) 



Figure 2: A sample of two topologically different Riemannian polyhedra hav- 
ing the same spectral properties. It is stronger than the (n— 1)— chainability 
and is aimed at avoiding topologically different Riemannian polyhedra which, 
however, can have the same spectral properties 

Let first 7 G QndAi be a boundary (n— I)— dimensional simplex with its 
(n — 1)— dimensional interior denoted by 7 mt . We introduce boundary normal 
coordinates in an relatively open subset U C Q as 



Here a(p) = d(p,dA4), and we assume that there is a unique q G 7 int with 
d(p, dM) = d(p, q), such that p lies on the normal geodesic to q, <^(r), <r„(0) = 



q, S v {d(p, DM)) = p and q v (Q, d(p, DM)) G ft int . If s(q) = (a 1 , . . . , s (n_1) ) are 



some (local) coordinates on 7, e.g. baricentric coordinates, then s(p) = s(q). 

Let now 7 C fi_ H fi+ be an (n — 1)— interface between n— simpleces Q- 
and Q + . Let U± be relatively open subsets of Q± with the nearest point 
on dfl± lying on 7 such that (U_f > |7) = (U+f^). Denote by (s, cr±) the 
boundary normal coordinates in U±, where s = (s 1 , . . . , s^ n-1 ^) are some 
local coordinates on 7, e.g. baricentric coordinates with respect to f2_ or 
Q + . We introduce the interface coordinates (s,a) on U_ IJU+: 



p^ (s(p),<t(p)), pet/. 
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Then the metric element in these coordinates takes the form, 

(dlf = {da) 2 + {g ± ) af 3{s, a) ds a ds p . (3) 
Throughout the paper we assume that the following condition takes place: 

Condition 2.2 For any interface 7 and any point q = (s, 0) on 7, the metric 
tensor g a p has a jump singularity at q. 

2.3 Laplace operator 

Let H l (M) be the Sobolev space of functions u G L 2 (J\4) such that Ui = 
u\Q i G and, for any interface 7 between f2j and flj, 

J -y I 'Y * 

The inner product on H l (M) determines the closed non-negative Dirichlet 
form, 

1 

B[u, v] = ^(ui, Vi) H i {ni) , u, v G H l {M). 

i=l 

By the standard technique of the theory of quadratic forms, the form D 
determines a self-adjoint operator in L 2 (Ai), namely, the Laplace operator 
with Neumann boundary condition, A. The domain T>(A) is defined by 

T?(A) = {ue H\M) : B[u,v] = (f,v) L , (M) for some / G L 2 (M)}, (4) 

where v G ^(M) is arbitrary Analysing condition (j4]), we see that u G 
T>(A) if u G H 2 (Qi), i = 1, . . . , I and, on any interface 7 C Q- f] Q+, 

?i_| 7 = n + | 7 , [ v /^ra (T M_]| 7 = [ v /#7<9 CT w+]| 7 . (5) 

where g±(s) = det[ {g±) a p{s, 0)]. 

As, due to the finiteness of Ai, the embedding of into L 2 (A4) is 

compact, the spectrum of A is pure discrete, 

= Ai < A 2 < . . . , A fc — > 00, 

with the corresponding basis of orthonormal eigenfunctions to be denoted 
by {tpk}kLi- Standard considerations, see e.g. [9] or [43J show that {<Pk}kLi 
distinguish points in A4 mt , i.e. for p 7^ q G .M int , there is k with (pk(p) 7^ 
<Pk(q)- Moreover, 
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Proposition 2.3 Let {y>k}kLi be an orthonormal basis of eigenf unctions 
of the Laplace operator A. Then {<Pk}kLi form local coordinates near any 
p G A4 mt , i.e. there are fci(p), . . . , k n {p) such that ((p^, • • • , ^Pk n ) form local 
coordinates near p. 

3 Gaussian Beams near interfaces 

3.1 Gaussian beams on smooth manifolds 

In this section we briefly recall some results on the non-stationary Gaussian 
beams on smooth manifolds. Their theory goes back to the pioneering works 
[5], [19], [ID]. In our exposition we follow mainly section 2.4 of [21]. Non- 
stationary Gaussian beams are some (formal) solutions of the wave equation 

U tt - AU = 0, (6) 

which are concentrated, at each moment of time t, near a point x(t). The 
point x(t) moves with a unit speed along a geodesic on a smooth Rieman- 
nian manifold (A/ - , h) with A being the Laplacian corresponding to (A/", h). 
Introducing a moving frame 

y{t) = x- x(t), 
a formal Gaussian beam has a form as a formal series 

U e (t, y) x M £ exp {-(ie)"^*, </)}$>(£, y)(ie) 1 . (7) 

l>0 

Here M £ = (ite)~% , < e < 1; 6 and u h I = 0, 1, are formal series in 
powers of y. They are usually represented as sums of homogeneous polyno- 
mials in y with coefficients depending on t, 

6 x ^6> m (i, y), ui = ^2ui m (t, y), 

m>l m>l 

6 m and ui m being homogeneous polynomials on y. The polynomials 9 m and 
M/ m are chosen so that, considered as formal series with respect to y and (ie), 

d 2 t U £ - AU £ = 0. (8) 

Note that " x " exactly means that the formal series ([7]) satisfies formally 
equation §8§. 

The most important properties of the non-stationary Gaussian beams are: 
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(a) 9i(t,y) = (£(t),y(t)) = £j(t)t/ J (t), where £j(t) is the unit covector corre- 

sponding to the geodesic x(t); 

(b) 02{t,y) = (H(t)y,y), where H{t) is a symmetric matrix, satisfying 

Im(H(t)y,y) > C(T)\y\ 2 , for -T < t < T. 

Remark 3.1 From now on throughout this paper we use the following nota- 
tions C (or, Ciit)) is a generic constant, C > 0, independent of e; fj,(L) is 
defined for sufficiently large positive integers L such that n{L) — > oo when 
L — > oo. 

Conditions (a) and (b) imply that U decays exponentially outside an e 1 ^ 2 - 
neighborhood of x(t). It is important to note that, starting from a formal 
Gaussian beam U £ we can construct a family of solutions to the wave equation 
(P), which "looks like" U e . To this end, we start with a finite series 

L 

U L e = M e exp{-(isr 1 e\t,y)}J2^(t,y)x(d 2 (x,x(t))8- 5 / (i ), (9) 

i=o 

L L 

6 L = S 9 m ] u\ = £ Ui m , 
1=0 1=0 

where x( s ) is a smooth cut-off function equal to 1 near s = 0. Then 

\\d 2 U £ L - AU £ L \\ c ^ HMx[ _ T , T]) < C L (T) e-« L \ 
By standard hyperbolic estimates there exists a solution U £ to (jSJ) such that 
IKE^-^Ho^xhT^i) < C L (T)e-^ L \ (10) 
Moreover, if we generate a wave inside M by a boundary source 

U £ \dAfx[-T : T] = fe(t,s). (11) 

Let f E (t, s), s G dj\f, t G [— T, T] be given by a formal expansion 
f e {t, s) x expi-iie)- 1 ^, s)}Y, 

l>0 

where 

§(*, S) X -t + £ a S a + < H((S, t), (S, t)) > 0m(t, S)l <? a/3 (0)£a£/3 < 1, 

m>2 
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ui{t,s) x y~]ui m (t,s), (12) 

m>0 

with 9 m , u\ m being homogeneous polynomials of degree m with respect to 
(t,s). Assume that ImH > 0. Then there is a unique formal Gaussian 
beam U £ satisfying (jSJ), and the boundary condition (JIT]) for — to < t < t 
with some to > depending only on geometry of (Af,dAf,h). Moreover, 
the corresponding geodesic x(t) starts, at t = 0, from the point s = 0, into 
the (co)diversion (£ Q ,£ n ) with £ n = [g^Ca^] 1 ^ 2 geodesic starting at s = 0. 
We will refer to this result saying that we can guarantee a non-stationary 
Gaussian beam propagating transversally to dAf by a proper choice of a 
boundary source (for these and other results on non-stationary Gaussian 
beams see e.g. [20], [2T]). 

3.2 Gaussian beams at interfaces 

In this section we consider reflection and transmission of the non-stationary 
Gaussian beams from and through an (n — 1)— dimensional interface 7 be- 
tween two n-simpleces f2_ and Q + . As our constructions will be of a lo- 
cal nature we can, without loss of generality, restrict them to the reflec- 
tion/transmission of the Gaussian beams from a smooth interface inside a 
smooth manifold. These questions, for the incidence angle less then criti- 
cal, were considered in detail in [25] and [39] . in the latter restricted to the 
isotropic media. Assuming that x(0) G 7 and introducing the interface nor- 
mal coordinates (s,a) with s = corresponding to x(0), we have, for the 
incident Gaussian beam, U™ 

00 

£/f(M)Uo x M e exp{(ie)- 1 e in (t jS )}5];(ie) , «r(t, a ) (13) 

1=0 

00 

[^gZd a Ui n (t,s)] Uo x M £ ex P {(z £ )- 1 §-(i, S )}^(z£)^(t !S ). 

i=-i 

Here O m , w| n , Sj n are sums of homogeneous polynomials with respect to (t, s) 
with 

B m (t, s) = -t + CV* + H((s, t), (s, t)) + ... ; 

Ao = C<, C>o, g a(3 CtT + (Zn) 2 = i, 
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with (Ca, CT) being the (co)direction of the Gaussian beam at t — 0. When 
(0)£^ n £g n < 1, it is possible to construct two formal non-stationary Gaus- 
sian beams t/J and U £ r in Q_ and Q + , correspondingly such that 



ur + u r £ , 



tr 



in fL 
in fij 



(14) 



satisfies the wave equation ([8]) and transmission conditions ([5]). To this 
end we use the technique briefly described in section 13.11 which reduces 
the problem to finding boundary conditions of form (fTTI) . (fl2l) for U r e and 
i7* r at a = 0. In turn, this is possible utilizing transmission condition ([5]) if 
(0)£^ n £g n < 1. Summarizing considerations of [25], we obtain the following 
result 

Lemma 3.2 Lei £7* n be a formal non- stationary Gaussian beam which hits 
the interface a = at s = 0, t = with its (co) direction (CaiCp) satisfying 
9+^(Q)€a€a n < 1- 27ien t/iere are two formal Gaussian beams U r e in and 
U £ r in Q + such that the total wave U £ satisfies the transmission condition 
(TJJ). The ( co ) directions (£o)£n) and (Ca^n) of the geodesies, corresponding 
to Ul and U £ r satisfy at t = 0, s = 0, a = the equation (SnelVs Law): 



tr rtr pin tr pin ttr r-i 



o/3/ 



\ £tr£tr\l/2 



The main amplitude coefficients u r 00 and u 
at t = 0, s = 0, a = by 



tr 
0,0 



of U r £ and Uf' are related to 



u 



0.0 



tr 
0,0 



'9- 



- cm 



7/*" ■ 



0,0 



trW 



(15) 



where g± = det 



We note that transmission condition for E/J n + f/J, and f7* r is under- 
stood in the formal sense. Namely, U T e and [7f r may be expressed in the form 
~g± d a U £ ,tr |(j = o having decomposition of form ( flBl) . Then 



(ED with u r / r \ a ._ 

(IBT) means that 



6' 



e r = e 



tr 



and 



it 



.tr 
7 ) 



sp + si 



Ui 
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as polynomial with respect to (t, s) 



Observe that the condition 



1 n-l 



¥> 9?{0) 



a,P=l 



n-l 



a,p=l 



implies that 



u r 00 ^ for almost all (CCD- 

When dealing with an incoming non- formal Gaussian beam (Q, (fTOj) . 
which we will denote by U™' L , similar to the above we find 6 r J ,i and [/* r,L 
by formulae (G| with 6 r,i , u{ L and Q tr,L , u\ r ' L instead of 9 m ' L , u] n,L , cor- 
respondingly. Clearly, they give use to approximate transmission conditions 

■"W. (17) 
Add to J7| r,i a function ^(t, s, a), 

L 



\\m n ' L + - U^ L )}\\c^,(-t ,t 0)) < C(T) e-> ,iL >. ( 16) 
\[(^gZdM n ' L + U r £ > L )) - V^d a Ul r > L ]\\ o » mbx( _ t0>t0)) < C(T) e 



k=0 



where ^k(s,t) are chosen so that 



U: 



a,L 



satisfies 



£ 1 £ 



[yfgId a (Av(U™> L + Ur> L ))] 



APTTtr,L\ . 



75 



for < p < [§]. Clearly (JTSD , (PU) imply that 



(18) 
(19) 



lcfW(n + x(-to,to)) 



< C(T) e 



Together with (1181) . ( 1191) . it follows from the wave equation (jHj) that we can 
modify and f/f ' L - 



r.L 



tr.Z, 



with 



wdr'm^M) <c L (t )e- 

\A^ L ^\\ LHM) <C L (t )e-^; 
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so ioi —to < t < to. Therefore, 

L _ jui n ' L + U^' L , a<0, 

satisfies the wave equation (jSJ) and coincide with {J l £ n,L for negative t, more 
precisely, for t < —Ce 5 ^ 12 . Here $^ is given by for a < and $* r>L for 
a > 0. 

Remark 3.3 When 7 zs (n — 1)— simplex in dA4, we can modify the previous 
construction to find the Gaussian beam reflected from 7. The part of Lemma 
EH related to U T e remains valid with formula for the main term u r Q taking 
the form 

u r 00 = -u™. (20) 

Thus, by the considerations similar to the above it is possible to find solutions 
to the wave equation fiEj) which satisfy the Dirichlet boundary condition and 
look like a Gaussian beam. 



4 First n— simplex 
4.1 

In this section we start proving the uniqueness Theorem 11.31 Recall that 
we are given diffeomorphic open subsets T C 7, T C 7, where 7_and 7 
are boundary (n — 1)— simpleces of n— simpleces Q C Ai, Q C M.. We 
assume, after proper unitary transformations in finite-dimensional spaces 
corresponding to eigenvalues of higher multiplicity, that 

\ k = A fc , (f k \r = <fk\f, fc = l,2, ... , (21) 

where x : F — > F is a diffeomorphism. Our goal in this section is to prove 
the following result 

Lemma 4.1 Let x : F — > F , F Q <s 7, r d 7 be a diffeomorphism satisfy- 
ing conditions of Definition \1.2\ Then there is an isometry X : Q — > f2, such 
that 

<Pk\n = ^*^fcln! fc = l,2,..., X|p = x. (22) 
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To prove this lemma, observe first that if To <e T, To d T and r > satisfy: 
exp an : T x [0,r) ->• M, exp^ : f Q x [0, t) ^ M 

are regular, i.e. map into Q and f2 correspondingly and provide an diffeo- 
morphism between T x [0, r) and exp 9n (r x [0, r)) C Q and T x [0, r) and 
exp 9n (r x [0,r)) C Q and, in addition, 

d(T ,M\n), d(f ,M\n) > r, 
then the diffeomorphism 

X = k x I : r x (0, r) -> f x (0, r) 

satisfies 

X Vfe|f x(o,r) = yfc|r x(o,T)- 

Moreover, X is actually an isometry between T x (0, r) and T x (0, r) 
considered as domains in (Q mt , g) and (Q mt , g), correspondingly. The proof 
of this fact is identical to the smooth case and is given in Section 4.4 of [21]. 

Assume now that Ui, u 2 C £l mt and Wi, w 2 C £l mt are open subsets with 
Xi : Ui — > Ui being isometries satisfying 

^»*Wfc|si = ¥>*L> i = 1, 2, k = 1,2, ... . (23) 

By Proposition 1, -STil^no^ = in<^2 which makes it possible to extend Xj 
into an isometry X : u)\ U a;2 — > ui\ |J £2 which also satisfies (1231) . 

Consider the family of all pairs of open subsets u C Q and a; C Q which 
are isometric to each other and satisfy ( 1231) . Clearly, this family is partially 
ordered by induction, by the above we can consider its maximal element 
which we denote by (ft m , fl m ) with the corresponding isometry denoted by 
X m . We want to show that Vt m = Vt mt , Vt m = f2 mt . 

To proceed, recall the following result from |21j . which is proven for 
smooth manifolds but remains valid for Riemannian polyhedra under the 
conditions formulated below. 

Theorem 4.2 1. Let S C Q (or S C Q) be a smooth subdomain such that 
{ { Pk{p)}k > ( or {'?kip)} ( k') are known for p G S (orforpE S). Assume 
that for r > 

expgs : dS x (0, r) —> Q \ S, exp d § : dS x (0, r) — > Q \ S 
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are regular. If, in addition, r < {d(S,A4 \fi), d(S,M then these 

data determine uniquely fk\s T y Pk\s ■> where 

S T = S\Jexp ds (dSx(0,r)), S T = S\J exp gS (dS x (0, r)). 

2. Let u, u be solutions of the initial boundary value problem, 

utt- Aw = F G C™(S xi + ); « tt -Au = FeC oo (5xl + ), (24) 

u\ t=0 = f G C °°(5); £| t=0 = / G C °°(^), 
wtltM) = G C °°(S); utl^o = G C °°(S). 

Then these data determine u, u on S T x R + , S x IR + , correspondingly. 
In particular, if S and S are isometric, with isometry X satisfying (H3j], 
there is an extended isometry X T , 

with 

¥k\s T = K&k\s T , u f \s T xR+ = X;u f \s rXR+ , (25) 
when f = X* f, (p = X* (f), F = X* F . 

4.2 

Based on Theorem 14. 21 we will finish the proof of Lemma I4~T1 Assume, that a 
maximal element (fi m , £l m ) ^ (£l mt , Q mt ), where without loss of generality 
we can take Q m ^ Q rnt . Therefore, there is a point p G Cl(Q m ) H Q mt . 
Observe, that as X m : Q m —>■ Q m is an isometry, it may be extended to the 
mapping X m : Cl{VL m ) — > CZ(fi m ). Consider the following possible scenarios: 

1. p = X m (p) G Vt mt . Denote by 5 = mm(d(p, dQ), d(p, dtt)) and by 
p = min(i(f2, g), i(Q, g)), where i(N, h) stands for the injectivity radius 
of the normal coordinates of Riemannian manifold (N, h) . Let 5q = 
±min(<5,p) and p G Q m , p = X m (p ) G fi m satisfy d(p ,p) < S . Let 
< cr < S satisfies B(p ,a) C Q m , B(p ,a) = X m (B(p ,a)) C Q m2 
where B(p, r) is a closed ball of radius r centered at p. Taking S, S 
in Theorem 14.21 to be B(p ,a), B(p ,o~), we see that conditions of this 
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Theorem are satisfied for r = 25 with S a = B(p ,a + 2<5 ), .S CT = 
S(po 5 0" + 2#o)- Therefore, X m can be extended to Q m U B(p , a + 2c> ) 
containing p, which contradicts the definition of Q m . 

2. p = X m (p) G 7 m *, where 7 is some (n — 1)— subsimplex of f2. Let now 
S = min(d(p, dQ)), d(p,d(Q U f^i)), where Qi is another n— simplex 
adjacent to 7 (if 7 G <9.M we take d(M\Q) rather than d(QUQi)). As 
earlier, let So — \ min(<5, p, i^(p)), where i^(p) is the radius of injectivity 
of the interface normal coordinates related to p, (or to the boundary 
normal coordinates of p G dM). Introduce po, po = X m (p ) as in the 
case 1 and take balls B(p ,a), B{p ,a), < a < 5 similar to the case 
1. Consider now the non-stationary Gaussian beams Uf 1 ' 1 on M. which 
start at t — at po, in direction close to the normal direction from p to 
7, i.e. with the initial co-vector £ close to (0,0, These Gaussian 

beams reflect from 7 and return to B(p , a) approximately at the time 
t = 2<i(po,7) — cr- Thus, for £ close to (0, ...,0, 1), the total Gaussian 
beam = U^ L + U l £ n > L in satisfies 



when q G B(p , cr), t G [2d(p , 7) - a, 2 d(p , 7) + a). 

On the other hand, the corresponding Gaussian beam [/* n ' L in Q moves 
from B(p ,a) and has no reflected part for t G (0, 2rf(p 0l 7) + cr)), so 
that 



When e is sufficiently small and L is sufficiently large, (1261) . fl2"71) con- 
tradict the second equation (|2^|) . of Theorem 14.21 As (fk\B(p ,a) an d 
(fk\B(p ,a) are known we can evaluate the Fourier coefficients Uk{t) = 
u k (t) of U^(t) and U^(t) and, next, 

U^q,t) = Y,Uk{t)Vk{q)\ U^(q,t) = £ u fc (t)&(g) , 

for g G B(p ,a), q G B(p Q ,a). 

3. Let now p G .M™ -2 . Using e.g. baricentric coordinates in Q, we see that 
there are C > 1 and 5 > such that if d(p~o, dQ) < 5, then there is a 
curve x : [0, 1] — ► f2 m *, x(0) = po 5 ^(1) = §0 such that: 



max|L/" e L (c/,t)| > C, 



(26) 



max |[7" £ L (<?,*) I < C^ (L) , for g G B(p ,cr). 



(27) 
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(i) d(x(t),dQ) > C^d^dQ); 

(ii) \x[0, 1]| < C d(pa, dft), where \x[a, b]\ is the arclength of x(t) be- 
tween x(a) and x(b); 

(iii) d(q ,dQ) = d(po,dQ). However, we can assume that there is a 
unique nearest point q G <9f2 to q and, in addition, 

(iv) go £ 7 mi ) where 7 is an (n — 1)— interface between Q and some Qi 
(or 7 C dM); 

(v) Interface (boundary) normal coordinates centered at §0 are regular 
in 4 d(qo, g)-vicinity of g - 

Returning to the consideration of the case p G M. n ~ 2 , let p G fi m 
satisfy d(p ,p) < min(5, ^5), where 5 is defined as in the case 1. With 
Po = X m (p ), let x(t) be a curve in Vt mt described earlier. Then the 
diffeomorphism X" 1 can be extended onto some open neighborhood of 
x(t). Indeed, by the construction of the step 1, we can move recurrently 
along x(t), using the balls of radius ^d{po, dQ). By the maximality of 
(Q m ,Q m ), q G Vt m with q = X" 1 ^) G Q m satisfying 

d(q ,dn)>Ad(q ,q). (28) 

Inequality (128]) makes it possible to use the same considerations as in 
the case 2, proving that the case p G M. n ~ 2 is possible. 

We finish the section by the following 

Corollary 4.3 Let X : Q —>■ Q, where Q, Q are n-simpleces of M. and Ai, 
satisfy / f^j) . If 7 = X(j) are (n — l)—simpleces of Q and Q correspond- 
ingly. Then 7 is an interface if and only if 7 is an interface. 

Proof Let / = / e L , (p — 0^ be the initial data for Gaussian beam JJ l £ n,L 
starting in Q towards 7 and (f^, (fi^) = X* (/^ , (p £ ) be the initial data of 
the Gaussian beam Uf 1 ' 1 . Comparing U£' and U^ L = X* U^ L and using 
formulae (fToT) and (I20I) . we obtain the desired result. □ 
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5 Global Isometry 
5.1 

In this section we show that isometry on X : Q —>■ Q which satisfies fl22|) can 
be extended to a global isometry X : AA — > AA satisfying 

Mp) = X*Mp), k = l,2,..., P eM mt ,p = X(p)eM mt , (29) 
d(p,q) = d(p,q), p = X(p), q = X(q), p,qE AA. 

We start with the following result which is a partial generalization of Theorem 
EES 

Lemma 5.1 Let 7 be the interface between Q- and Q + in AA. and^y = X^(j) 
is the interface between and Q + in AA. Assume that X_ : — > is 
an isometry satisfying |Hj]. Then X can be extended to an isometry X : 
n_Ufi+^fl_U £1+ which satisfies (TUP on Q™* U fi£*. 

Proof Let D, D be smooth subdomains in Vt™ 1 and Vt™ 1 with the "upper" 
part of the boundary parallel to 7, 7, i.e. the parts of dD, dD are given 
by x = (s, — o"o), p = (s, —<To), s G T d 7, s G r d 7. Assume, without 
loss of generality that d(D,d(Q- U d(D,d(Q- U > 10cr , a < 

|(min(zr, if)), where ir, if are the injectivity radii of interface coordinates 
related to T, T, correspondingly. We want to show that equation fl22l) implies 
that 

(Pk(s,°~) = V?A:0s,<t), for s G T, s = X_| 7 (s), — 2cr < o < 2<r . (30) 

We have the following rather straightforward generalization of Tataru's ap- 
proximate controllability: 

Clv iM ){u f (to),fe C™(D x (0,t ))} = L 2 (U t0 (D)), < t < 4a , (31) 

with a similar identity for D. Here U to (D) is a to - neighborhood of D in AA 
and it^ is a solution to the initial boundary value problem 

u{ t -Au f = f, u f \ t<0 = 0, u f \ dMxR+ = 0. 

As usual this result follows immediately from observability, see e.g. Section 
2.5 of [21]. To formulate the desired observability, let 

v tt - Av = 0, in M x (-t , ^o), «|aMx(-to,to) = °> 

v\Dx(-to,t ) = °, 
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where < to < 4<t . Then v = in double cone 

K(D,t ) = {(p,t) eMx (-t , to) : d(p,D) <t - \t\}. 

To use this we first observe that for any q with d(p, D) < 4cr and any 5 > 
there is a piece-wise smooth curve x : [0, 1] — > U Q + with x(0) = p G 
.D, x(l) = g and \x[0, 1]| < d(q,D) + 5. Moreover, x(t) may be chosen to 
cross 7 transversally. Then we can continue v by along this curve x(t) so 
that for x(s) 

v = in K x (-(t - |ar[0,s]|, t - |x(0,s)|)), 

where V s is a small vicinity of x(s). Indeed, this is obvious for pieces of the 
path x(s) lying inside either Q 1 ™ 1 or To cross 7 we just observe that, if 
v, d a v = on £ x (— t, £ ), £ (e 7, when approaching 7 from Q + , then by (8), 

v, d a v = 0, on E x (—£,£), 

when approaching 7 from f2_ , and, therefore, v can be continued by further 
along x(t) into Q + and vice versa from Q + into Clearly, to use (jSJ) we 
should assume t> to be sufficiently regular, however, by smoothing v with 
respect to time t, we can extend it to non-smooth solutions, see e.g. [22] . 

This implies that v(q,t) = for \t\ < t — d(q,D) — 5. As 5 > is 
arbitrary, we obtain that u = in K(D,to). 

Identity (jHTj) makes it possible, starting from {A^, ^Id}^, and from 
{A fc , ^15)^=1, such that A fc = A fc and (fk(s,cr) = tp k (s,a), for -2<x < 
o < —o~o, J = (X-( , y)(s)), s G T to construct (fk(s,a) and (pk(s,a) for 
2<x < a < 2a, s = (X_(7))(s), s G T, see [20], Chapter 4.4 of [21]. In par- 
ticular, the construction in [20], [21] imply that ipk(X-.( i y)(s),a) = ipk(s,a). 
Observe now that T x [0, 2<To], T x [0, 2o~ ] from a relatively open subdomains 
in fl + , Q + , respectively with X + : (s, a) — > ((X_| 7 )(s), a), < a < 2a, being 
a diffeomorphism satisfying (I2"5j) . 

Mimicking the proof of Lemma 14.11 in section I4.2[ we extend X + to be 
a diffeomorphism, X + : Q + — ► + , satisfying (122]) . As, by construction, 
X + | 7 = X_| 7 , X defined as X_ on f2_ and X + on Q + is a desired isometry 
between f2_ U f2+ and f2_ U 

5.2 Identification of n-simpleces 

We have proven that, for any chain of n— simpleces Q\ = f^m, fii( m ) m 
Ai which are pairwise adjoint, i.e. there is an (n — l)-interface 7^ be- 
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(M,g) (M,g) 

Figure 3: Chains of chambers, connecting points 

tween fl^k) and Q^k+i), k = l,...,m — 1, there is a chain of n-simpleces 

Qi = fij(i), fii( m ) in .M, such that fij(fc) and fij(fc) are diffeomorphic with 
diffeomorphism X^) satisfying (1221) (strictly speaking Xj(fc) may depend on 
a chain from Q\ to ^i(fe)), see Fig[3j Let us show that 

(a) For any n— simplex Qj C M. there is an n— simplex Qj C M. which 

is diffeomorphic to Qj with diffeomorphism Xj satisfying f)22p and, 
likewise, for any Qj C M. there is Qj C M. with described properties; 

(b) If f2j and f2j are diffeomorphic with diffeomorphisms Xj and X'- satis- 

fying d22]) then Xj = X' j ; 

(c) If ilj is diffeomorphic to fjj and fl. by Xj and X^- satisfying f[2"2"j) . then 

Let now fij(i), ^j(i) be f2, f2 are the ri— simpleces with 7, 7 used in Lemma 
14. II being their boundary (n — 1)— subsimpleces. Then, by Lemma [4.11 there 
is Xi : fij(i) — > %(i) satisfying (1221) . Denote by 71 the (n — 1)— interface 
between fij(i) and Qj ( 2 ) - By Corollary 14. 3[ there is f2j(2) adjacent to Qj(i) 
with the (n — 1)— interface 71 = fij-m H %(2) sucn that 71 = Xi(7i). By 
Lemma EU1 there is an isometry X2 : ^(2) ~ ^ ^3(2) with Xi| 7l = X2| 7l . 
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Continuing this process, we obtain an isometry X m : £lj( m ) = ty, ^j(m) := 
which satisfies (1221) . 

To prove (b), let p G Qf \ and p G Xj(jp), p = X'^p) G tif*. As X j7 X] 
satisfy (I2"2"j) . ^(p) = <Pk(p), k = 1,2, .... By Proposition 12. 3[ p = p . As 
p G fij 1 ' is arbitrary, Xj = X- on Using the fact that {<Pk}kLi distinguish 

points in J\A mt , see Proposition 12.31 we prove property (c). 
Based on properties (a)-(c) we show the following result. 

Lemma 5.2 Let LBSD for the Riemannian polyhedron (A4, g) and (Ai, g) 
be equivalent. Then 

1. For any Qi in M. there is a unique diffeomorphism X, : f2* n * — > Qf 11 , 
which satisfies l[2~2]). 

2. 7 is an (n — 1) — interface between Q_ and fl + if and only if the corre- 
sponding 7 is an {n — 1) —interface between f2_ and Q + . In this case 

X_| 7 = -X+| 7 , 

where X_ is the closures of the described above diffeomorphisms X T on 

3. The diffeomorphisms UX can be uniquely extended to an isometry 

X : M reg -> M reg . 

Here X is an isometry of A4 reg and A4 reg considered as metric space 
with the distance function given in Definition \2.1\ 

Proof By (a)-(c) it remains to prove the part of this lemma dealing with 
(n — 1)— simpleces. Let 7 be an interface between Q_ and Q + with Q_ and Q + 
being the corresponding n— simpleces in Ai. Crossing 7 we move from fL_ to 
U in M. so that, by the previous constructions, (l mt is a diffeomorphic to Cl + 
with diffeomorphism satisfying (1221) . But also f2+ and f2+ are diffeomorphic 
with diffeomorphism satisfying (|22|) . Thus f2 = Vl + and the diffeomorphisms 
X T : fl!?* — > are uniquely extendable to an isometry X_ !+ : Vt mt U 
Qint y ^ _^ Qint u Qint y w h ere 7y — X(7). As the distance functions in 

A4 re9 , M. reg employ only curves being in M\M n ~ 2 , M\M n ~ 2 , the above 
results concerning the local isometries yield the desired global isometry. 
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Theorem 11.31 immediately follows from Lemma 15.21 taking into account 
the Definition 12.11 of an admissible Riemannian polyhedron. Note that this 
definition implies that the topology of a Riemannian polyhedron considered 
as a metric space is the same as the topology of the underlying polyhedron. 

6 Further generalizations and estimates 
6.1 

Condition 12.21 that the metric tensor does have a jump singularity at every 
point of any (n — 1)— dimensional interface, i.e. Condition 12.21 may be too 
restrictive. In this section we relax it a bit. namely, we assume that it 

• either the metric tensor g does have a jump singularity at any point 
p g ynt^ f or an ^ n _ i)_i n terface 7, between and Q + , 

• or g is smooth across 7 at any p £ 7*™*. 

The latter condition means, that in the interface coordinates related to 7, see 
section 1.2, the metric tensor g a p{s, a) is smooth at o = 0. In this case we call 

7 an artificial interface and would like to treat f2_ Uf2 + together, introducing 
f^*U7 m *Uf^. Further removing artificial interfaces and taking into account 
that Riemannian polyhedron consists of a finite number of n— simpleces, it 
is natural to introduce the following object. 

Definition 6.1 A chamber Q of a Riemannian polyhedron (Ai, g) is the 

maximal union of open n— simpleces together with open artificial interfaces 
adjacent to them. 

It is natural to treat each chamber Q as an open Riemannian submanifold of 
M. with a piece- wise smooth boundary. However, due to [21], there may be 
topological obstruction to that. Namely, it may happen that, for p £ AA n ~ 2 , 
there is no open neighborhood of p (in Ai) which is homeomorphic to R n . In 
this connection we introduce 

Definition 6.2 An admissible Riemannian polyhedron (A4,g) is called weakly 
admissible if any of its open chamber Q mt is an open n— dimensional Rieman- 
nian manifold with piece-wise smooth boundary and, if (n — l)—subsimplex 
7 of M is not artificial, they either 7 £ dAA or the metric g has a jump 
singularity across 7. 
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Then, similar to Theorem 11.31 we can prove 

Theorem 6.3 Let (Ai, g) and (Ai, g) be weakly admissible Riemannian 
polyhedra. Let local boundary spectral data corresponding to (Ai, g) and 
(Ai, g) are equivalent. Then there is an isometry X : Ai — > Ai such that 
for any open chamber f^ n * C Ai, X\ Q int is a diffeomorphism onto an open 

chamber Qf 11 C Ai, which satisfy equation [2^1 . 

Local boundary spectral data in this case consist of a "smooth" open subset 
T C dAi, such that there is an open neighborhood U C At, of To, which 
is diffeomorphic to a half-ball {x G B(0, r) : x n > 0} with Tq being an open 
subset of B(0,r) n {x n = 0}, and eigenpairs {Xk, d v <pk\r }kLi- 

6.2 

Similar to the smooth case, the uniqueness Theorem 11.31 remains valid when, 
instead of the boundary spectral data, we have a local non-stationary Dirichlet- 
to-Neumann maps, 

AF : C°°(T x (0,T)) - C°°(r x (0,T)), 

where C°°(r x (0, T)) consists of smooth functions which are qual to near 
t — 0, is then defined as 

A r : / 9 v u f \ To x{o,T), 
where u* is the solution to the inverse boundary value problem 

U>tt-Au f = 0, U f \ t<0 = 0, U f \dMx(0,T) = /• 

Condition 12.21 may be replaced by a condition that the metric tensor g 
is not smooth across any (n — 1)— interface 7 at any point p G 7 m *. This 
condition means that, for any p G 7 m * and any interface conditions with 
s = corresponding to p, the metric tensor g a p(0, o) or some its derivatives 
d^g a p(0,a) has a jump at a = 0. Then Theorem (11.31) remains valid as, 
for regions of 7 where g a p{s,a) is continuous across a = together with 
its first (k — 1)— derivatives with respect to a, but d^.g a p(0,a), does have a 
jump across a = 0, the reflected Gaussian beam U r e is of order k. Namely, 
in representation for U^(t,y) the first non-zero ui(t,y) is Uk(t,y). This, 
however, does not alter considerations of section [3]l5] with the only difference 
that in section 15.11 the estimate for the reflected Gaussian beam in (1261) is 
changed to | U r £ ' L (q, t) | > Ce k . 
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6.3 Open problems 

(i) This paper deals with uniqueness in the inverse boundary spectral prob- 

lem for Riemannian polyhedron. It does not provide an algorithm to its 
reconstruction. In particular, trying to apply the technique of Section 
4 [21], when approaching M. n ~ 2 the size of a step of the reconstruc- 
tion procedure tends to zero. Therefore, it is impossible to reconstruct 
(Ai, g) by a finite number of steps. 

(ii) Even with generalization described in sections 5.1, 5.2, the class of Rie- 

mannian Polyhedron considered in this paper does not cover an impor- 
tant case when the metric g is smooth across some part of the interface 
7 but is not smooth across other part of 7, i.e. we have "holes" in 
interfaces. We intend to study such cases in the forthcoming paper. 

(iii) Another important open question is the one of stability which could 
relate observation error with , in addition to curvature, injectivity radii, 
etc., [4j, size of n— simpleces, value of jumps, etc. 
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